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Abstract. Approximate analytical closed energy formulas for semirelativistic 
Hamiltonians of the form cr^/p 2 + m 2 + V(r) are obtained within the framework 
of the auxiliary field method. This method, which is equivalent to the envelope 
theory, has been recently proposed as a powerful tool to get approximate analytical 
solutions of the Schrodinger equation. Various shapes for the potential V(r) are 
investigated: power-law, funnel, square root, and Yukawa. A comparison with 
the exact results is discussed in detail. 
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1. Introduction 

Since the beginning of quantum mechanics, there has been a considerable amount 
of works devoted to the computation of analytical solutions of the Schrodinger 
equation, especially in bound states problems. Apart from its intrinsic mathematical 
interest, finding analytical formulas is always useful in physics, for example to 
obtain informations about the dependence of observables on the parameters of the 
Hamiltonian. This can be a great advantage when one tries to fit the parameters of a 
model to some experimental data. 

Many methods allowing to find approximate analytical solutions of the 
Schrodinger equation for bound state problems have been extensively used in the 
literature: WKB method, variational method, etc. (many solved problems can be 
found in p] for example). Recently, we have proposed a new method to compute the 
bound states of a given Hamiltonian (2j [3]. This method, called the auxiliary field 
method (AFM), is based on auxiliary fields - also known as einbein fields [4] - and 
can lead to approximate closed analytical results in many cases: One can quote central 
potentials of the form ar x + br v for special values of A and rj [3] , exponential potentials 
of type —ar x e _/3r [5], and square root potential [6[[7]. 

In standard Schrodinger Hamiltonians, the kinetic term is of the form p 2 /2m. 
However, there exists another class of Hamiltonians in which the kinetic term is 
relativistic, typically 

H = ay/p 2 + m 2 + V(r), (1) 

with a a positive number. Such Hamiltonians are usually called "spinless Salpeter 
Hamiltonians" . They are actually a well-defined approximation of the Bethe-Salpeter 
formalism applied to the description of bound states within quantum field theory [5] . 
The underlying assumptions is that the constituent particles interact instantaneously 
and propagate like free particles [9]. Moreover, the spin effects are not taken into 
account in this framework. For this reason one often speaks of "semirelativistic 
Hamiltonians" to characterize the spinless Salpeter Hamiltonians. 

Because of the square root in the kinetic energy, it is no mystery that finding 
analytical eigenvalues for Hamiltonians like |T]) is even more difficult than in the 
nonrelativistic, Schrodinger, case. Such a task can however be fulfilled by using 
envelope theory (ET) [10j . which is a powerful method to get approximate analytical 
solutions of eigenequations, for various potentials [HI [TSJ [13l HU [15] . We have shown 
in a previous work that AFM and ET are, to a large extent, two equivalent approaches 
[IB] , That is why it is reasonable to suppose that, using the AFM, we will be able 
to get analytical formulas for the eigenvalues of spinless Salpeter-type Hamiltonians. 
That is the topic the present work is devoted to. 

Our paper is organized as follows. The general principles of the AFM are recalled 
in section [2[ and its application to spinless Salpeter Hamiltonians is discussed in 
section [3] The method is applied to analytically solve several cases of interest in 
physics: Power-law potentials in sections [4] and (5[ square root potential in section [6l 
funnel potential in section [71 and Yukawa potential in section [8[ By comparing our 
formulas to numerical results, we then show in section [9] that their accuracy can be 
considerably improved by performing some minor formal modifications. Finally we 
sum up our results in section 1101 Since the formulas we obtain are valid for one- 
particle systems or for two particles with the same mass, we give in | Appendix A| some 
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results in the unequal mass case, while some polynomial equations that are of interest 
in this work are presented in | Appendix B| 

2. The auxiliary field method 

2.1. General technique 

The auxiliary field method has been first presented in [5] and extended in [31 [5] . In 
order for the present paper to be self-contained, we recall here the main points of the 
AFM and refer the reader to the aforementioned references for more details. 

The aim is to find an analytical approximate solution of the eigenequation 
H\%) = E\$), with 

H = T(p 2 ) + V(r). (2) 

As usual, we work in natural units h = c = 1. The AFM consists in the following 
procedure. We first assume that Ha = T(p 2 ) + v P(r), where v is a real parameter, 
admits bound states with an analytical spectrum. The eigenequation 

H A \*{v))=E A {v)\*{v)) (3) 

is thus analytically solvable, as well as the one associated with the Hamiltonian 

H{v) = T(p 2 )+V{r,v) with V(r,v) = vP{r) + V{I{v))-vP(I{v)). (4) 

The function I(x) is defined by 

I{x)=K-\x) and K(x) = ^-, (5) 

P (X) 

the prime denoting the derivative with respect to x. Let us remark that V{r, v) is of 
the form C\P{r) + C2 where C\ and C2 are constants. These numbers are determined 
in order that V(r) approximates at best the potential V(r). Using © we have 

E(v) = E A (u)+V{I(u))-vP(I(y)). (6) 

The approximate eigenvalues and eigenstates are eventually given by E(i/q) and 
I^A^o)) respectively, with vq such that the total energy ([6]) is extremal, i.e. satisfying 

d„E(v)\ v=Vo = 0. (7) 

The value of Vq depends on the quantum numbers of the state considered. Once 
vq is known, the constants C% and Ci can be computed. It is shown in |16j that 
V(ro,Vo) = V(r ) and V'(r , Vo) — V'{ r o), where r is such that Vq = K(ro). 

The above procedure can be justified as follows. Let us assume that, instead 
of considering V(r, v) in which v is a constant, we consider V(r, v) in which v is an 
arbitrary function. A proper elimination of this auxiliary field by the following 



constraint ScV(r, C>) 



leads to the solution vq = K(r). The original 



Hamiltonian ((2|) can then be recovered since V(r, vq) = V(r). The essence of the 
AFM is then to replace this function vq by the constant v$. In [2], it is shown that 
vo ~ (^a(^o)I ^0 \^a{vo)}- The AFM can consequently be regarded as a "mean field 
approximation" with respect to a particular auxiliary field which is introduced to 
simplify the calculations. The main technical problem of the AFM is the determination 
of an analytical solution for ([5|) and 0. Such a task has already been fulfilled for 
nonrelativistic Hamiltonians, either with central potentials of the form ar x + br TI [3] 
or with the exponential and Yukawa potentials [5] . 
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An estimation of the error on an exact eigenvalue E can be obtained. In [5] , it is 
shown that 

E{v ) -E> V(r Q ) - (y A (vo)\V(r)\y A (v )). (8) 

The r.h.s. of this equation is the difference between the value of potential V computed 
at the "average point" rg and the average of this potential for the trial state. As we 
are mainly interested in the computation of eigenvalues, no evaluation of the quality 
of the approximate eigenstate I^aC^o)) nas been performed. This will be studied in a 
subsequent paper. 

An important result concerning the AFM has to be pointed out: This technique 
is equivalent to the envelope theory [10j . which is an other method to compute 
approximate analytical energy formulas [16 . As a consequence of this equivalence, 
the energy formulas which are found within both frameworks are identical. Since 
many results concerning semirelativistic Hamiltonians have already been obtained by 
resorting to envelope theory [T21 E] one could wonder if it is still relevant to perform 
the present study. We claim that the answer is positive, mainly because in [12l 114) . 
the energy formulas are all written under the form E = min w [/(f)]. The minimization 
stage ([?[)■ which is generally non trivial, is thus not explicitly performed with the 
consequence that no closed analytical formula is given; we achieve this stage in the 
present work in order to have completely explicit formulas. Moreover, we consider 
in the following two semirelativistic Hamiltonians that have not been studied up to 
now: the ones with square root and Yukawa potentials, both of interest in theoretical 
physics. A last argument is that we also propose an improvement of the energy 
formulas thanks to a comparison with the exact numerical data as it is done in [21 [21 [5]. 
So, we are finally led to more accurate analytical formulas. 



2.2. Useful results 

We have shown in [2] that the eigenvalues of 

H x (a) = El+asgn(X)r\ (9) 
v 

where sgn(A) = A/|A| is the sign of A and where a is a positive parameter, can be 
written, by using the AFM, under the form 



e A (a, v, N) 



2 + A 


aX 


2 

A + 2 




A 


T 




(v) 



(10) 



N depends on the radial n and orbital I quantum numbers. An important result, 
proved in [3J, is that the functional form of (fit))) does not depend on the potential 
P(r), provided that this potential is of power-law form. Only N actually depends on 
the particular form of P(r). With the choice P(r) — r 2 , N = 2n + I + 3/2. Using the 
results from [TB], it is then possible to show that (Ti"T)|) gives an upper bound of the 
exact result. With the choice P(r) = — 1/r, iV = n + / + 1. A lower bound is then 
yielded by (fT0|) . Since the form of N is an artifact of the AFM, it can be modified to 
better fit the exact results. In particular, we have shown that, for values of A £ [—1, 2], 
a good form for N is given by [2] 

N x = b(X)n + l + c(X), (11) 

with 

, m _ (4^-18)A + (18-2c) (7c-36)A+(36-5c) 
1 ' (3w-15)A+(21-3w)' 1 ' (6^-32)A+(40-6w)' 1 1 
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and w = i/3ir. Such expressions for the functions b(X) and c(A) lead to energy formulas 
as precise as 10~ 3 for the most interesting (the lowest) values of the quantum numbers 
n and I [2]. But the variational character of (fT0|) cannot then be guaranteed. Notice 
that the energy formula (|10p together with (jTTj) and (TT^J) gives the exact result for the 
harmonic (N2 = In + I + 3/2) and Coulomb (iV_i = n + I + 1) potentials and also 
reproduces the exact asymptotic behavior of the linear case (Ni = -Z=n + l + ) [2] . 
We also found in [5J an accurate analytical energy formula for the Hamiltonian 
P 2 b 

H = — + ar . (13) 

v r 

This linear-plus-Coulomb potential, also called the funnel or Cornell potential, is 
of great importance in hadronic physics. It corresponds to a linear confinement 
associated with a short range Coulomb contribution coming from one-gluon exchange 
processes |T7j , and has been shown to emerge as the effective energy between a static 
quark-antiquark pair in lattice QCD [18j . Defining 



one can express the energy spectrum of (| 13[) as [3] 

ef =V3^(-^_--l-) , (15) 



Fl(Y) F + (T), 

with F + (Y) given by (TR2]) . With the choice N = N 2 (N-i), (H3) gives upper (lower) 
bounds. A better accuracy can be achieved with the choice N = b((3) + 1 + c((3) , where 
/3 4 = v 2 b 3 /(27a). The explicit form of b{(3) and c(/3) can be found in [3]. But again 
no information on the variational character can then be obtained. 

As we have analytical and accurate energy formulas for any power-law potential 
and for the funnel potential at our disposal, we will assume in this work that the 
choices P\(r) — sgn(A) r x or Pf(r) = ar — b/r are good starting points to apply the 
AFM. 



2.3. Alternative method 

There exists another way to introduce auxiliary fields that has already been presented 
in [19) ; we recall and extend it in this section. This alternative method is formally 
simpler but less general than the AFM presented in section^ Nevertheless, it can lead 
to analytical results in many relevant cases and will appear to be of crucial importance 
to find analytical eigenvalues of semirelativistic Hamiltonians. 

Let us suppose that the Hamiltonian of our problem is the sum of Q different 
contributions 

Q 

H =Y,K h (16) 
t=i 

where Ki are some operators. The eigenvalues and eigenstates of Hq are denoted Eq 
and \ipo) respectively. We can introduce k auxiliary fields (k < Q), denoted as (f>i, to 
obtain a new Hamiltonian Hk as follows 

H k (<f> m ) = + Pi<f>i) + E 

i=l V ^ / ]=k+l 
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with real numbers a,, /3j such that c^/?,; = 1/4. Since the idea is to replace an 
operator by its square, this method can be expected relevant only for definite positive 
operators. In this case, the quantities a*, ft, and 4>i must also be definite positive. 
The auxiliary fields 4> m should be rigorously eliminated as operators thanks to the 
constraint 5d, i Hk(<j>m)\j,—£- = 0- One then obtains 



<f>i = y/at/PiKi, (18) 

and Hk((f>j) — H$. Consequently, Hamiltonians (]16p and (|17p are equivalent up to a 
proper elimination of the auxiliary fields as operators. 

We now consider that tfo are variational parameters in order to simplify the 
calculations and denote Ek(<fi m ) and \ipk{4>m)) the eigenvalues and eigenstates of 
Hk(4>m)- Then the set of coupled equations d ( j >i Ek((p m )\^ =r p — have to 
be solved and the final eigenenergies and eigenstates read Ek-o = Ek{4> m -.o) and 
|"0fc;o) = \i>k{<l>m\o)) respectively. An example where this procedure can be useful is 
the linear potential: No analytical solution can be found in general for the Schrodinger 
equation with V{r) — ar, but the replacement of V(r) by aV 2 (r)/(f) + @<j) [sec (fl~7|) ] 
only involves a harmonic oscillator and an analytical solution is then possible to find. 
Another application which will be used below is the elimination of the square root 
operator in (pj [20] . 

The Hellmann-Feynman theorem [5T] states that 

dfcEkfym) = (V>k(0m)| d^H k ((j) m ) \lpk(<f>m)) ■ (19) 
The computation of this equation for the fields <pi-o gives 

4>l,0 = ^ (V'fciO I Kf |^ ; o) = (A-,o I # l^o) ■ (20) 

Pi 

It is then possible to obtain an analytical form for the approximate energy 

k Q 

E k , = A&;0 + £ (V>fe;o| #j MkiO) 

i=l J=fc+1 



= Ev^°i^ 2 i^+ E fci^iv^o)- (2i) 
i=i j=fc+i 

The result is then independent of ai and fy, which can be chosen at best convenience. 
At first sight, this equation could be seen as a triviality, but the problem is to determine 
\tpk-,o}- This is precisely the purpose of the AFM. In the rest of this section, we 
will denote (ipk-fll O \ipk-fi) by (O). Let us remark that ([20)) confirms the mean-field 
interpretation of the AFM. 

This eigenenergy Ek-o can be rewritten as 
k 



Using the trivial inequality 



2 



Ki - ^-<pi ) > 0, (23) 



2q 

together with aj/3j = 1/4, it is easy to show that 



Jr<<aA+/3 ( &. (24) 
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So, we have H k (cf) m ) > H. From [55] and [55], we can conclude that (H ) > E for 
any state. It follows that 

E k - fl >E +S with 5 = (\[(KT) - ( K i )) • (25) 

It is easy to see that 

< ^ (k - sfiKf)^ = 2^(kT) (sJiKT, (Ki?) , (26) 

which shows that S > 0. This implies that the energy spectrum obtained by 
introducing auxiliary fields, as described by the above procedure, is an upper bound 
of the exact energy spectrum. A lower bound on the difference E^.q — E is then given 
by 5. 

3. Semirelativistic Hamiltonians 

The rest of this work is devoted to find analytical approximate energy formulas for 
semirelativistic Hamiltonians of the form 

H = o-^/p 2 + m 2 + V(r), (27) 

where the potential is assumed to be central, and where the real parameter a can take 
all positive values. Actually, a — 1 or 2 are used for most practical problems. If we 
deal with a one-body problem in a central potential or if we consider a two-particle 
system with a particle of mass m and a particle of mass M 3> m (so that the dynamic 
of the massive body can be forgotten) , we must take a = 1. On the other hand, in the 
two-body case, we restrict our study to systems where both constituents have the same 
mass m and, in this case, a = 2. Comments about the unequal mass case are given in 
|Appcndix A| Although these two situations are the most usual, we find convenient to 
let a as a free real parameter; some examples of such a useful prescription will be given 
later. The relativistic nature of this last Hamiltonian comes from the fact that the 
kinetic energy has no longer a nonrelativistic form, i.e. p 2 /2m, but a relativistic one. 
The operator \/p 2 + m 2 is actually the relativistic kinetic energy of a free spinless 
particle of mass m. The neglect of spin effects is one of the reasons why we use the 
word "semirelativistic" to characterize Hamiltonian (|27p. In the literature, such a 
Hamiltonian is usually referred to as a spinless Salpeter Hamiltonian. 

3.1. Auxiliary field formalism 

Because of the square root appearing in (|27[l . finding exact analytical eigenvalues 
seems to be hopeless, excepted in very particular cases (see section I5"T2"|) . However it 
is possible to resort to the alternative AFM presented in section 12.31 to rewrite the 
spinless Salpeter Hamiltonian under an apparently nonrelativistic form. Let us set 
K\ = a ^p 2 + m 2 and K 2 = V(r) in (TBJ). Then, we can introduce one auxiliary field, 
denoted as <j>\ — v, with a.\ = 1/er 2 and j3i = o- 2 /4:, and rewrite (| 1 T[) in this particular 
case as 

2 2 2 

H(u) = — + ^-u + ^-+V(r). (28) 
v 4 v 
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The auxiliary field v can be eliminated as an operator thanks to the constraint 
S v H(u)\ u=j> = 0, leading to 

0= -vV + m 2 . (29) 

(T 

Obviously, one has then H(v) — H, that is the equivalence between Hamiltonians (f27|) 
and ()28j) . But, if v is seen as a real number, (|28|) is simply a nonrelativistic Hamiltonian 
with kinetic term p 2 /is, for which analytical solutions can be hoped. To our knowledge, 
this technique has firstly been used in [53] but without explicit reference to the AFM. 
If e(v) are the eigenvalues of h(v) = p 2 jv + V(r), one can write the eigenvalues of 
H{v) as 

E(v) = ^ + + e(u), (30) 
and the final spectrum is given by E(i/q), where vq is such that 

2 9 

9uE(u)\ u=Uo = => — + e> ) = — . (31) 

The prime denotes a derivation with respect to v. Provided that e[y) is analytically 
known, the only difficulty in this procedure is the analytical resolution of this last 
equation, especially when to ^ 0. Reporting the value of vq given by (|3T|) in the 
expression of the energy ([30|) allows to write 

2 

E{v ) = ^v Q + {ve{v))'\ v=VQ . (32) 

With the notations used in this section, E is an exact solution of the Hamiltonian 
H (|2"T]) and E(v) is an exact solution of the Hamiltonian H{v) (f2"5|) . From the results 
of section [2~3l we know that E(v) > E for any state, and that E(vq) > E, since it 
is also true for the particular value vq of the auxiliary field. Let E*(v) be a solution 
(often approximate) of the Hamiltonian (|28|) obtained by the AFM. If E*(i>) > E(v), 
then we have E*(vq) > E, and an upper bound is available for each eigenstate of the 
Hamiltonian (I2T)) . If E*{u) < E{v) or if it is not possible to locate E*(u) with respect 
to E(v), no information can be obtained for the variational character of E*(i/q). The 
value of isq is given by 



v Q = -^{p 2 + m 2 ), (33) 

G 

where the mean value is computed with an eigenstate of EL(vq). An estimation of the 
difference between the exact and approximate eigenenergies is given by (|25[) . which 
reads in this case 



6 = o-^(p 2 + to 2 ) - a ^ Vp 2 + m2 J ■ (34) 
At the limit of low mass, (|2T|) can be written 

9 

fT 'IT I i 

H^H UI + —== with H m = ay/p^ + V{r). (35) 

The contribution A (to) of the mass to to an eigenenergy of the ultrarelativistic 
Hamiltonian _ff ur appears as a small contribution that can be computed as a 
perturbation. The mean value being taken with a eigenfunction of H UT , we have 

... / am 2 \ am 2 m 2 . „. 

A(m) = ( — = > fa — == = — , 36 

where Po is the auxiliary field obtained with the AFM applied to the Hamiltonian 
H m . The validity of this approximation will be checked in section [7731 
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3.2. Scaling laws 

Let us consider the spinless Salpeter equation 

(a vV + to 2 + G V(a r) ~ E(m, G, a, cr)) *(r ) = 0, (37) 
where G and a are two real numbers. Defining p — ar, ir = p/a, one has 

f V^ 2 + {m/af + — V(p) - —E(m, G, a, cr) ) *(p/o) = 0. (38) 
\ acr acr / 

This is just the spinless Salpeter equation 

(y/^+m 12 + G' V(p) - E'{m', G', 1, 1)) <p(p) = 0, (39) 

with 

m'=™, G' = -^, £V,G',l,l)^ (m ' G ' a ' g) . (40) 
a acr acr 

We are thus led to the following scaling law for the energy spectrum 

E(m,G,a,a) =aaE[ — , — ,1,1) . (41) 
\ a aa J 

It means that both one parameter of the potential and the number of constituents in 
the system can always be set equal to 1 in order to simplify the computations without 
loss of generality. The general energy formula will then be recovered thanks to (14"T]) . 

It is worth mentioning that in the particular case of a power-law potential, i.e. 
V(r) = Gsgn(A)r A , the same computations as above can be performed and lead to a 
more constraining scaling law that reads 

E x (m, G, cr) = (cr A G)^r Ex ( X = (g) ^ m, 1, l) . (42) 

This means that the whole power-law problem can be studied by only considering the 
eigenvalue E\(x, 1, 1) associated to the two-parameter Hamiltonian 



H x ( X ) = v 7 ^ 2 + X 2 + sgn(A)/ . (43) 



4. Power-law potentials 

4-.1. General case 

We apply in this section the AFM to find approximate analytical energy formulas for 
the spinless Salpeter Hamiltonian ([27]) with the potential V\(r) — asgn(A)r A where 
a is positive. We consider that A > —2, as needed in the nonrelativistic case to have 
bound states. 

First, we can get rid of the square root appearing in the kinetic energy by 
introducing one auxiliary field v as it is done in (|28|). We obtain 

2 2 2 

H\(v) = — — h ^—v + — + asgn(A)r A . (44) 
v 4 v 

If v is assumed to be a real number, this last Hamiltonian is simply a nonrelativistic 
one with a power-law potential, for which analytical energy formulas were obtained in 
[2] and recalled in section l2~2l Using the result (fT0|) one is led to 

2 

N^. (45) 



_ , , m 2 a 2 A + 2 _ 

V 4 A 



with A 



x 



a A 
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At this stage, setting N = N\ seems a good choice (see (fTTj) 1 ). but it is not necessarily 
the best one in the present case. Indeed, N\ has been fitted on the energy levels of 
the exact nonrelativistic spectrum, while we deal here with a relativistic kinematics. 
Thus, we can expect that the best expression for N is different from N\. Inspired by 
(fTTj) . one could try N = b(m, X)n + I + c(m, A) (see section |9j). 

The auxiliary field v can now be eliminated by minimizing (|45|l with respect to 
v, the minimum being obtained for the value v = vq. Setting 

A + 2 

v Q = x 2 , (46) 

we obtain 

d,E x {v)\ v=Vo = => ^4+ 2 - A\ x - m 2 = 0. (47) 
Assuming that we know the solution xq from (|47|) , two convenient forms can be found 
for the physical values E\ = E\ (x 2 ) of the eigenenergies (|45l) 
a Am 2 + (A + 1) A\ x 

Ex = \ / 2 i a ' ^> 

A sjm z + A\ Xq 

or 

In order to obtain analytical energy formulas, one should be able to analytically 
solve (|47p . Let us now examine the cases for which this is possible. First, we set 
A + 2 = p/q and Xq = x^ q , where p and q are relatively prime. Then, ()47|) becomes 

^XP-A x X q Q -m 2 = 0. (50) 

A polynomial possesses analytical roots if its order is less or equal to 4; therefore 
all the solvable potentials should verify the conditions 1 < p, q < 4. An exhaustive 
research of all the solvable potentials leads to the following values for the power A: 

A = -2, -I -|, -|, -1, -5, -!,-§, -I, 0,1, 2. (51) 

Among these allowed values, three are of particular interest: A = — 1 and 2 because 
the Coulomb problem and the harmonic oscillator play a central role in theoretical 
physics, and A = 1 since a linearly rising potential is generally considered to be a 
relevant approximation of the confining potential in QCD |17j . These three cases are 
explicitly solved in the following sections. 

Let us briefly discuss the existence of physical solutions, as function of the 
parameters. Starting from (|4"T| and (|4"5f . it is possible to prove the following properties: 

• For A > 0, there exists always a physical solution with E\ > 0; 

• For — 1 < A < there exists a physical solution with E\ > as long as the 
parameter a is less than a critical value a c (A) given by 

-M-'U) w G£f ■ 

Let us remark that a c (X) depends on m, except a c {— 1) = a N. 

• For — 2 < A < — 1 either there is no solution for (|47|) . or there exist two solutions 
but in the latter case these solutions are not compatible with the nonrelativistic 
expressions when m tends toward infinity. In both cases the corresponding 
solutions are not physical so that these types of potential must be discarded 
from our study. 
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4-2. Harmonic oscillator 

If A = 2, (|4"7} becomes a quartic equation of the form 

2 

^-Xq - A 2 x - m 2 = 0, with A 2 = y/aN. (53) 

Defining 

f2A 2 \ 1/3 m 2 ( 16a \ 2/3 , s 

(|53p can be rewritten as 

4X 4 -8X- 3Y 2 = 0, (55) 

which is precisely of the form (|B.7[) . Following (|B.8|) . the solution of (|55|) is given 
by X(Y 2 ) = G-(Y 2 ), and, after a rearrangement of (|48|) . the energy spectrum reads 
(using (|55p under one of the equivalent forms) 

E 2 =amJ— , 56 



3/2 K 



2 



r 2 VG-(y 2 ) 2G^(r 2 ) 



(57) 



2cr m /„ g , 1 

G 2 _(Y 2 ) + 7r ^ 7T . (58) 



V3y V ~ v ' g_(f 2 ) 

With N — N 2 these formulas yield an upper bound of the energy, since this setting 
gives the exact solution for the corresponding nonrelativistic Hamiltonian. 
In the nonrelativistic limit (Y — > oo), the last equations reduce to 



E 2 &o-m+\ —N, (59) 
V m 

as expected for a nonrelativistic harmonic oscillator (see (jTUJ) for A = 2 and v — 2m /a). 
For large value of m, the choice N = N 2 is clearly optimal. For small values of m, 
another choice could give better results. 

4-3. Linear potential 

The resolution of the case A = 1 is rather similar to the one of the harmonic oscillator. 
Defining 

AiV /2 v , „ 3 3 / 2 m 2 c 



Xo = 2 {^) X and Yx ^^TT' m 

(|47p simply becomes a cubic equation of the form (|B.1|) . that is 

X 3 - 3X - 2Yi = 0. (61) 
Notice that A\ = (aN/2) 2/3 so that 
SVSam 2 

2 aN { ' 
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Following (|B.2|) , the solution of ([CT]) is given by X(Y{) = F_(Yi), and, after a 
rearrangement of (|48|l . the energy spectrum reads (using (|61[) under one of the 
equivalent forms) 

/T" iy+3JL(yo 

L 1 =am\ — ; (63) 

V*i y/3F_(Y 1 ) + 2Y 1 V 7 

= <rmj - . ^— ) (64) 



(3 + F?(yO). (65) 



In the nonrelativistic limit (Y — > 00), the last equations reduce to 



3 /aa 2 N 2 



1/3 

fii*<rm+5(— I , (66) 
2 \ m / 

as expected for a nonrelativistic Hamiltonian with a a linear potential (see (|10p for 
A = 1 and v = 2m/ a). For large value of m, the choice N — Ni is a good one. For 
small values of m, another choice could give better results. 

4-4- The Coulomb problem 

Equation (|4"7|) is considerably simplified when A = — 1; the value of xq can directly be 
extracted in this case and reads 



Am 2 . , , / a \ 2 

-0 = ^3^, with A_x = (-) . (67) 



The energy spectrum is finally given by = my/ a 2 — 4A_i, or 



^_ 1= amy 1-^5. (68) 

It is obvious from this last equation that bound states exist only when a < a N, in 
agreement with (|52[) . The most stringent upper bound for a is actually found for 
the ground state, that is a < a N\ n=l=Q . With N = AT_ X , |[B5]) yields an upper 
bound of the energy since this setting gives the exact solution for the corresponding 
nonrelativistic Hamiltonian. 

As the coupling constant a is dimensionless in this case, the only mass scale of the 
Hamiltonian is m. So the nonrelativistic limit cannot be obtained by setting m — ► 00, 
as usual. It is well known that the mean speed of the particle is independent of m and 
proportional to a. So, the nonrelativistic limit is achieved for a — > 0. When a <C 1, 
(|68j) reduces to 

„ a 2 m 

E ^ am ~^ (69) 
which corresponds to the rest energy plus the Coulomb binding energy, as expected 
(see (fTU)) for A = —1 and v = 2m /a). For large values of m, the choice N = N-i is 
clearly optimal. For small values of m, another choice could give better results. 

Several works have been devoted to the spinless Salpeter equation with a Coulomb 
potential [231 HS1 123 HZ] > and it is interesting to compare our results with previously 
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found ones. First of all, it has been shown in [53] that the energy spectrum of the 
semirelativistic Coulomb problem is unbounded from below if 

a> a c = — . (70) 

Moreover, a lower bound on the ground state energy E g is [24] 



E g >am^l--. (71) 

But, since we used the alternative AFM to get rid of the square root in the kinetic 
term, (|68[) is an upper bound of the exact energy. This is confirmed by the results of 
|26j in which formula (|68p was already given. 

A last result of interest is the analytical determination of the ground state energy 
performed in [27] . It is found that 

Eg\ a = ac < am X 0.4842564... (72) 

Our formula (l68l) leads in this case to 



£_i| 0=0o =<rmJl- — — 2 ■ (73) 

V n Ar ln=i=0 

This last expression is equal to 0.77m if N — N-x, but agrees with (172"]) if N\ n=l=0 m 
0.73 is taken (see 19.21 for more information). 

5. Ultrarelativistic limit 

In the ultrarelativistic limit, that is to say m = 0, the various elimination equations 
for the auxiliary field and the energy formulas become simpler. Before going further, 
we stress that the results presented in this section are only valid for m = 0. 

5.1. Power-law potentials 

Equation (|47[) can be written in the ultrarelativistic limit as 

x (^a$ +1 -A>) =0. (74) 
The nontrivial solution of this equation is 

and the final energy spectrum is given by 

E^(a,a, N) = |aA|^r ( a N)^ . (76) 

A 

Let us stress that, in this special case, an analytical expression is obtained whatever 
the power A of the potential, in contrast to the particular values (|5T]) resulting from 
the general case. One can check that (l56l) and ([63]) reduces to ([76]) when m — ► 0. 

For the particular case A = 1, this last formula becomes E™' = 2 V a aN, and 
one has consequently (i?" 1 ) 2 oc N. Such a linear behavior between the squared 
eigenenergies and the quantum numbers of the different states (with N = N\ for 
instance) is a well-known property of the spinless Salpeter Hamiltonian with a linear 
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potential and massless constituents. Such a Hamiltonian is one of the simplest ways 
to describe a light meson in a potential approach of QCD, and it is experimentally 
checked that the squared masses of the light mesons mainly grow linearly with their 
angular momentum (Regge trajectories). See for example [17] for a discussion of that 
point. 

The quantity EY is physically a mass. It appears that EY < when — 1 < A < 
and ET^i — 0. No bound state of massless constituent particles can be found in these 
cases. Although ([76)) is positive for — 2 < A < —1, this range of values for A has been 
proved to be unphysical. 

By comparison of the nonrelativistic result (|10[) and the ultrarelativistic one ([76]) , 
it can be checked that 

e 2X (a 1 ^,arh/a,VN) =E?(a,<r,N). (77) 

There is a kind of duality between the nonrelativistic and ultrarelativistic energy 
formulas for power-law potentials, at least for the approximate relations given by the 
AFM. 



5.2. A special case: the Harmonic oscillator 
The energy spectrum of the Hamiltonian 

ffho = 2 v /p I + ar 2 (78) 
can be analytically computed for I = and reads [33] 

E ho = - (4a) 1 / 3 a„, (79) 

where a n < are the zeros of the regular Airy function [n = 0, 1, . . . ). They can be 
found for example in [29] . where the following WKB approximation of these zeros is 
also given 



3tt / 3 X x 2/3 



2 V 4 



(80) 



This approximation is precise up to 8%. With the WKB expression (|80|) . an accurate 
approximation of the energy levels ([79]) reads 

\ 2 /3 




E u ,, :-! ( v „ ( ^n+^Jj . (81) 

This last expression can be compared to the AFM result ([75)) for A = a = 2, i.e. 

E^ = 3(Va~N) V3 . (82) 

With ([75)1 as a starting point, a natural choice for N is N 2 , but it is clear that both 
.Eho and E™ are identical provided that 

N\ l=0 = ^n + ^. (83) 

These values for the coefficient of n and the zero point energy are exactly those that 
optimize the AFM energy formula in the case of a nonrelativistic kinetic energy with 
linear potential (see (O). The explanation of this fact is that the Fourier transform of 
Hamiltonian (|78p is a nonrelativistic Hamiltonian with linear potential. Consequently, 
it is not surprising that we find the expression l|83p for N\ l=0 . This is one more 
indication that N should be different in the semirelativistic and nonrelativistic cases. 
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6. Square root potential 

The square root potential V(r) — v a 2 r 2 + b 2 plays an important role in the theoretical 
description of hybrid mesons [30]. In the flux tube picture indeed, an excitation of 
the gluonic held gives rise to an excited flux tube potential of such a square root 
form. That is why we find interesting to show that analytical energy formulas for this 
potential can be easily obtained with the AFM. 



6.1. Nonrelativistic kinematics 

2 

The Fourier transform of Hamiltonian H = ^IVa 2 r 2 + b 2 , denoted as F(H), reads 



T{E) = o-^p 2 + M 2 + K r 2 , (84) 

with 

4a\ 1/3 .. b 



,2 



ma 

M = —, n = <7—. (85) 
a 8 



This last Hamiltonian is nothing else than a spinless Salpeter one (j2"T)l with a harmonic 
potential. Provided that the proper substitution rules are taken into account, the 
energy spectrum in this case is the same than the one computed in section [4721 We 
have thus 

E- = bSl ^+ 4G -^) ( 86 ) 

sq V y 2 v /8G_(r 2 ) + 3r 2 v ' 

^/Jf^ + ST^W) (87) 



Y 2 \G.(Y 2 ) 2G 2 _(Y 2 



2b f Gl(Y 2 1 



w 2 V g_(f 2 ; 



b 2 /32m^ /3 



with 

F2 =3^J ■ (89) 

This potential has been studied in more detail in [SJ [7j , and we recovered here the 
results presented in that works. It can be checked that when b — > 0, these last 
equations reduce to the energy spectrum of a nonrelativistic Hamiltonian with linear 
potential, that is (fT0|) with A = 1. 



6.2. Ultrarelativistic limit 

The Fourier transform of H = a yfp 2 + V a 2 r 2 + b 2 reads 



T(H) = ay/p 2 + M 2 + ar, (90) 

with M = b/a. This last Hamiltonian appears to be a spinless Salpeter one (|2"T|) with 
a linear potential. Provided that the proper substitution rules are taken into account, 
the energy spectrum in this case is the same than the one computed in section 14.31 
We have thus 

2~ yj+3F-(Yi) 

Yi v /3F_(yi) + 2y 1 



E™ = b x l-^===^= (91) 
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with 




(92) 
(93) 



la eriv 

When 6^0, these last equations logically reduce to (|76|) with A = 1. 
7. Funnel potential 

As we stressed in section 12.21 the funnel potential V(r) = ar — b/r, with a and b 
both positive, is of crucial importance in hadronic physics. It is thus an interesting 
potential to be studied with the AFM. The energy spectrum of ([28]) with the funnel 
potential is analytically known thanks to (I15|) . and is given by 



with 



Y = - and V = 3« 2 ,&. (96) 

V V 

Note that an alternative expression, often more convenient, for f|95[) is obtained thanks 
to the following relation 

FfiY) Ejyj = 2 { F+{Y) Ejy)) ' (9?) 

It is more useful to express the energy as a function of Y rather than as a function 
of v, since the analytical minimization with respect to Y looks formally simpler than 
the same task with v. After some algebra, the condition dyE(Y)\ Y — Yo — can be 
rewritten as 

V3ab m 2 a 2 U 



0. (98) 



3Fl(Y ) U AY 2 
The property F+(Y) + 3F + (Y) - 2Y = has to be used 

7.1. Massless particles 

In the ultrarelativistic limit (m = 0), (|98|) becomes 
F+(Y ) = _2b_ 
Y 3aN 

whose solution is 



D, (99) 



To obtain this result, the definition relation of F+(Y) is used with the special value 
Y = Yq. It is worth mentioning that one has the constraint 2— 3D > in order for Yq to 
be defined, or equivalently b < aN. This constraint is formally the same as for the pure 
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Coulomb potential and states that b cannot be arbitrarily large without forbidding 
bound states. This feature of the ultrarelativistic funnel case has numerically been 
observed in [3"T] . 

Finally, the energy formula reads 

£ f ur = Ef(Y ) = 2^/a(aN -b). (101) 

It is astonishing that such a complicated potential used in conjunction with a 
semirelativistic kinetic energy admits analytical approximate eigenenergies of such 
simple form. For 6 = 0, (| 1 1[) reduces to 2 V aaN, that is the expected expression from 
(I76|) in the case A = 1. For a = 0, E UI — > as shown in section [5j If the prescription 
N = N2 is chosen, then (|101j) is an upper bound of the exact result. As already 
mentioned, other choice for N could improve the accuracy of the formula but without 
no guarantee about the variational character of the result. 



7.2. General case 

At first sight, solving ([98]) . may seem hopeless. However, one can indeed obtain an 
analytical solution. The calculations are somewhat complicated, and we report here 
only the main steps and useful results. 

The first step consists in multiplying this equation by F + (Y ) 3 = 2Y — 3F + (Y ) 
(this is the result (|B.2[) applied to the value Y — Yq) in order to get an expression 
linear in the F function. One is able to obtain the following relationship 



where we introduced intermediate quantities 

2m 2m b 3 , 
aU 60N 1 V 3a 

A = aN 2 -bm 2 , (105) 

and where D has already be defined in . 

Instead of using the variable Yo, it is more convenient to work with the new 
variable Z defined as 

Z = cY 2 , (106) 

so that (|102[) can be recast into the simpler form 

, x 2Y (l-dZ ) , , 

with 

u 2 bra 2 

Let us remark the interesting relation c+3cd = D 2 which will be used in the following. 

Inserting p07[) into the equation defining the F + function (|B.2[) , we get rid of the 
F + function and we remain with a polynomial in Z . Explicitly, we end up with the 
equation in Zq to be solved 

4d 3 Z 3 + (D 2 - Y2d 2 )Z 2 + 6{D 2 + 2d)Z + 9D 2 - 4 = 0. (109) 
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Since the polynomial is of third order, its root Z can be expressed analytically. The 
same property holds for Y through (|106[) . for F + (Y ) through (|10T|) and lastly for the 
searched energy E through (|95p (the use of the alternative form (|57|) is interesting). 

Although we have obtained the analytical expression for this general case, the 
resulting form is quite complicated and not reported in this paper, but it can be 
obtained on request. We prefer to give the simpler case valid for low masses. 



7.3. Low mass expansion 

For m = 0, it is easy to show that u = d = 0, c = D 2 so that the solution of (|109p is 
Z = Z = (2 - ZD)/D leading to F + (Y ) = F = Vz = DY = DY, and E given by 
(jlOip . Let us define mo = N \J a/b. By low mass, we mean the condition m -C mo, 
or £ = (m/mo) 2 -C 1. In such a case, A — aN 2 (l — e) so that c = D 2 (l — e) and 
d = e(l — e) _1 /3. Keeping only the terms of first order in e, the equation to be solved 
(|109j) becomes 

L> 2 Z ( 2 + (6L> 2 +4e)Z + 9L> 2 -4 = 0. ( n 0) 
Since Z is the solution of this equation when e = 0, we search for a solution of 
the form Zq = Z + rj, with rj <C 1. Substituting this value of Zq in (jllpp and limiting 
ourselves to first order terms in e and 77, we find the solution 

z = z(i-Vi. an) 



D 



This allows to calculate 



Yo = Y (l + £^ e ) (112) 

and 

F + (y ) = ^(i-^). (us) 

Inserting these values in the expression ([93)) and (|97[) for the energy, we obtain 
the final result 



Here again, it is amazing that the eigenvalues of a so sophisticated Hamiltonian take 
such a simple form. The contribution A(m) defined by (|36[) can be computed with 
the value P derived from (J96J) , (JMl), (fTOO]) . One finds 



A W = wV— (115) 

Taking into account (jlOip . this is in perfect agreement with (|114| . 



8. Yukawa potential 

8.1. Energy spectrum 

To our knowledge, only a few analytical results exist concerning a relativistic Yukawa 
problem [321 133] . These results are moreover obtained by solving one-dimensional 
Hamiltonians, either Klein-Gordon or Dirac, but not a three-dimensional spinless 
Salpeter Hamiltonian. As we show in this section, the AFM allows to gain some 
insight about the latter case. 
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Let us start from (f2"T| with V(r) = —ac P r /r. With the change of variables 
x = (3r, q = p//3, one obtains 

n = — R = Vq~ T T^-g— , (116) 

ap x 

with 

ma _, 
.9--- H7 

We can moreover assume from section 14.41 that the eigenvalues of the Hamiltonian 
y/q 2 + 9 2 — g/x are exactly known and read 8^/1 — g 2 /N 2 . Consequently, we can 
apply the AFM with V(x) — —ge~ x /x and P(x) — —g/x. A quick calculation leads 
to 

v = <T I(v) [I{y) + 1] . (118) 
Such an equation can be inverted to obtain 51 

/(!/) = -1 - W_! (™) , (119) 

where W-± is the secondary branch of the Lambert function (or Omega function or 
product-log). The interested reader will find more informations about W-\{x) in [5]. 
Thanks to (JHJ) and (|119[) . the eigenenergies of (|116[) can be written as 



-Wi-^-sr^Tsr < 120 > 

The physical value v$ for v is the one ensuring the condition d v E y (v, g)\ v=v = 0; 
explicitly, one finds 



v c / x v q g * n z 

Notice that the property W'(z) = W{z)J [z(l + W(z))] is used to get this last relation. 
The insertion of this value vq into (jl20[) provides the approximate AFM energy 



X N 2 + X 2 gv Jl ^ 



Ey{y 0l g) = - — - y N2 -. (122) 

xgVo + N 2 ^Jl- 9 -^ 

The final energy spectrum in physical units is given by a E y (vo, g = ol/<j) with 
Ey(fQ,g) given by (|122[) and vq solution of p2ip . 

Unfortunately, (|121[) is transcendental so that v$ cannot be analytically extracted; 
this drawback is also present in the nonrelativistic version of this problem [5j. We 
nevertheless notice that in the limit (3 — > one has vq = 1 and (3 a E y (l,g) reduces to 
(|68|) as expected. 

It is possible to get an asymptotic expression for (|122p in the limit x — > oo. In 
order to simplify the notation, let us set e = N/x & small quantity and p = g/N. At 
the limit e — > 0, the third term of (| 1 2 1 [) vanishes and it is easy to check that vq — > 1 
by lower values. Thus, it is natural to set Uq = 1 — 77, with r\ a small positive quantity. 
Then, 1 + W-i(-i>o/e) -> -V2rj and the third term of (fT2T|) tends to e>/l - p 2 / p 
(neglecting the term proportional to rje). With these considerations, (|121[) reduces to 



y/2rj = eyl — p 2 / p, which immediately provides us with the rj value and hence the 
approximate Vq value. Explicitly, we have 

iV 4 / a 2 
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One deduces the approximate expression for the energy, valid for large values of \ 



E V M,X, 9) « 9 + (x - Q V 1 - j^- (124) 

In physical units (see (jTT6j) and (|117ll ). this last formula reads E y = a(3E y (v^x — 
m/(3,g = a/a), namely 



8.2. Critical heights 

It is well known in nonrelativistic quantum mechanics that the Yukawa potential 
—ge~ x /x admits a finite number of bound states. Consequently, there exists critical 
heights, denoted as g n i. They are such that if g > g no i , a bound state exists with 
radial and orbital quantum numbers n = no and I = Iq respectively. These critical 
heights have been intensively studied in the nonrelativistic case [321 |3S], but their 
computation remains problematic with a spinless Salpeter Hamiltonian [36) . 

Actually, one has to deal with two critical heights in a semirelativistic problem. 
To understand this fact, let us start from g = 0, for which E y (v, %, 0) = \ logically, 
and increase its value. Bound states progressively appear; in particular the bound 
state with quantum numbers n and l appears for g x - n i > such that 

E v( u x;0,X,9x;nl) = X^ v xfi g x , n i = . (126) 

X 

If g > g X ;nl indeed, E y will be smaller than the rest mass term x- Such a negative 
binding energy shows the existence of a bound state. Replacing i/q by v x -o in (11211) 
allows to obtain an explicit expression for g x -,ni, that is 



N 2 cV /ir ? 

fc "' = T^Tf <127) 

Note that the property W-i(xe x ) = x, which is actually a definition for the Lambert 
function, has to be used. In the nonrelativistic limit, x ~* 00 an d X9xM ~ * eN 2 /2 as 
expected from [5], where the nonrelativistic Yukawa problem is solved with the AFM. 
The factor 1/(2%) is indeed a consequence of the scaling laws in the nonrelativistic 
case, while the dependence in N 2 has been checked by comparison to numerical data 
0. 

We now continue to increase g. Since (|122[) is actually a mass formula, one cannot 
have a value of g so high that E y (vo,Xi g) < 0. But the numerator of (| 1 22[) can never 
vanish since x, v, N and g are positive. The only trivial possibility is x = 0, but it 
leads to an unphysical state with E y — 0. However, there is an obvious constraint 
stating that gv$/N < 1. Since vq G [0, 1] from (|118[) . one has 

9 < N. (128) 

Thus the energy ()122|) decreases with increasing g until g = N is reached. Then vq = 1 
is a solution of (|121[) and one has the minimal value E y (l, x, N) = N. No bound state 
with E y < N can exist. 

There are in consequence several constraints ruling the existence of bound states in 
the Yukawa potential. In order to have a complete spectrum, one must satisfy x > N 
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from (|127p . and g x . nl < g < N\ n=l=Q . However, the inequality g x . ni < N\ n=l=0 
is not trivially satisfied. One has g x ^ n i/N — f(N/x), with f(y) an monotonous 
increasing function of y € [0, 1] such that /(0) = and /(l) = 1. Thus, the inequality 
N/x < / -1 (-Nj n= i=o /N) should be verified. Finally, the constraints, expressed in 
physical units, are 

N a , .. 

m>P f-HN^/NY ^<-<^U=o- (129) 

The existence of the ground state is guaranteed if m > /3 N\ n=l=0 and g x -oo < &/o~ < 



9. Improved formulas 

The AFM provides numerous analytical approximate formulas for eigenvalues of the 
semirelativistic Hamiltonian |T]). They all depend on a global quantum number N. 
The AFM cannot give clear indications about the exact form of N. In previous papers 
El [5l [7] , it is shown that the accuracy of the formula can be greatly improved by 
fitting the form of N with the exact results. In this section, the same procedure will 
be implemented for some Hamiltonians considered above. Details are not given here 
but can be found in [21 [3] . It is worth mentioning that very accurate eigenvalues can 
be obtained numerically with the Lagrange mesh method [37j . It is very accurate and 
easy to implement, even in the case of a relativistic kinematics. 



9.1. Ultrarelativistic power-law potential 

Let us first consider the following dimensionless Hamiltonian 

H = 2y/t? + x x (130) 

with A > 0. The approximate energy spectrum is given by (|76p with a = 1 and a = 2. 
With the choice N = N 2 = 2n + I + 3/2, upper bounds are obtained. As mentioned 
before, another choice for the n- and ^-dependences of N can greatly improve the 
results. By using the form 

N = 6(A) n + d(X) I + c(A), (131) 

we find smooth variations for coefficients 6, c and d for A €]0, 2] (I < 3 and n < 3). In 
particular, d(X) w 1 with relative variations less than 2%. Coefficients b(X) and c(A) 
can be fitted with various functions and similar agreement. Finally, we choose 

3.00 A + 3.67 2.69 A + 8.69 

v ' A + 3.40 ' v ' A + 8.27 v ' v ' 

Agreement with exact results is very good but the variational character of the 
approximation is no longer guaranteed. With the choice (|132p . the maximal relative 
error for I < 3 and n < 3 and for A £ [0.1,2] is located between 0.3 and 1.1%. With 
the choice N — N2, the corresponding error is located between 4.5 and 12.7%. Results 
for A = 1 are presented in table [1] 

In section 15.21 it is shown that iV|._ must be given by (|8"3"]l for a quadratic 
potential. Formulas (|132[) give b(2) = 1.79 close to 7r/v3 » 1.81 and c(2) = 1.37 close 
tO 7T v^/4 » 1.36. This is also in agreement with formulas (71) in [2| which predict, 
in the case of a nonrelativistic Hamiltonian with a linear potential, 6(1) = 1.79 and 
c(l) = 1.38. 
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Table 1. Eigenvalues e(n, I) of the Hamiltonian H130H with A = 1, for some 
sets (n,l). First line: value from numerical integration; second line: approximate 
result with N defined by (1 1 3 1 I t and (1 1 32 I t : third line: approximate result with 



JV 


= 2n + I - 


f 3/2. 






I 


40,1) 


e (i,0 


e(2,0 


e(3,0 





3.f577 


4.7109 


5.8913 


6.8742 




3.1338 


4.6849 


5.8374 


6.7973 




3.4641 


5.2915 


6.6333 


7.7460 


1 


4.2248 


5.4575 


6.4837 


7.3767 




4.2215 


5.4725 


6.4866 


7.3623 




4.4721 


6.0000 


7.2111 


8.2462 


2 


5.0789 


6.1304 


7.0470 


7.8671 




5.0814 


6.1602 


7.0764 


7.8869 




5.2915 


6.6333 


7.7460 


8.7178 


3 


5.8108 


6.7425 


7.5775 


8.3387 




5.8156 


6.7785 


7.6207 


8.3787 




6.0000 


7.2111 


8.2462 


9.1652 



9.2. Relativistic Coulomb potential 

With dimensionless variables, the semirelativistic Coulomb Hamiltonian is written 



H = 2^+l-- (133) 

x 

with a c < a < 0. The approximate energy spectrum is given by (|68p with m = 1 and 
g = 2. With the choice N = N-i = n + I + 1, upper bounds are obtained and the 
results are exact in the limit a — > 0. As shown in the previous section, another choice 
for the n- and Z-dependences of N can greatly improve the results. By using the form 

N = 6(a) n + d(a) I + c(a) (134) 

and imposing b(0) = c(0) = d(0) — 1, we find smooth variations for coefficients b, c 
and d for a €]0, a c ] (I < 3 and n < 3). These coefficients can be fitted with various 
functions and similar agreement. Finally, we choose 

, . 1.03 o- 1.48 /N 1.07 a - 1.64 . . 0.96 a -1.56 , , 

6 <">= a -1.48 ' C{a) = a -1.64 ' = a -1.56 ' (135) 

Agreement with exact results is very good for a < 1.2 but the variational character 
of the approximation is no longer guaranteed. With the choice (|135p . the maximal 
relative error for I < 3 and n < 3 and for a £ [0.2, 1.2] is located between 0.005 and 
0.3%. With the choice N = N-%, the corresponding error is located between 0.004 
and 17.3%. To obtain a good accuracy in the domain a w a c — 4/tt w 1.273, special 
method must be used as the one presented in [57]. To recover the value obtained for 
the ground state in this paper, it is necessary to have c(a c ) = 0.73. Our formula gives 
c(a c ) = 0.77. Results for a — 1 are presented in tabled) 

One can see that N = N-i = n + 1 + 1 is a better choice for large values of n or I. 
This is can be understood as a kind of nonrelativistic behavior since the limits a — > 
and N — ► oo have similar effects. With / a generic name for the coefficients b, c and d, a 
more efficient form for these parameters should be /(a, n, I) with lim a ^o /( a i n i = lj 
lim ra _ >00 f(a, n, I) = 1 and limj^oo /(a, n, I) = 1. Such a refinement is beyond the 
scope of the present paper. 
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Table 2. Eigenvalues e(n,l) of the Hamiltonian JT33J with a = 1, for some sets 
(n,l). First line: value from numerical integration; second line: approximate 
result with N defined by 1 1134 \ and JT35J; third line: approximate result with 
N = n + l + l. 



1 


6(0,0 


6(1,0 


e(2,0 


6(3,0 





1.65817 


1.92184 


1.96739 


1.98231 




1.65982 


1.92356 


1.96680 


1.98151 




1.73205 


1.93649 


1.97203 


1.98431 


1 


1.93515 


1.97122 


1.98389 


1.98973 




1.93476 


1.97012 


1.98291 


1.98895 




1.93649 


1.97203 


1.98431 


1.98997 


2 


1.97187 


1.98416 


1.98987 


1.99297 




1.97296 


1.98416 


1.98961 


1.99266 




1.97203 


1.98431 


1.98997 


1.99304 


3 


1.98428 


1.98993 


1.99301 


1.99487 




1.98528 


1.99021 


1.99302 


1.99477 




1.98431 


1.98997 


1.99304 


1.99489 



9. 3. Ultrarelativistic funnel potential 

Written with dimensionless variables, the ultrarelativistic Hamiltonian with the funnel 
potential is given by 

H = lJtf +x-- (136) 

x 

with p > 0. The approximate energy spectrum is given by (|101[) with a — 2, a = 1 
and b replaced by P (to avoid confusion with the coefficient of n). This kind of 
Hamiltonian is often used in hadronic physics with typical values for P « 0.4. With 
the choice N = N2 = 2n + 1 + 3/2, upper bounds are obtained. As shown in previous 
sections, another choice for the n- and ^-dependences of TV can greatly improve the 
results. By using the form 

N = b{P) n + d(P) I + c{P) (137) 

we find smooth variations for coefficients b, c and d for P G [Q, 1] (I < 3 and n < 3); 
notice that P must be lower than the critical value 4/ir ss 1.27. These coefficients can 
be fitted with various functions and similar agreement. Finally, we choose 

1.88/3- 5.34 1.99/3- 4.40 ,,„. 0.76/3-2.46 

b ^= P-3M ' CW = P-3A9 ' dW = P-2M ' (138) 
Agreement with exact results is very good but the variational character of the 
approximation is no longer guaranteed. With the choice Q138p . the maximal relative 
error for I < 3 and n < 3 and for /3 e [0, 1] is located between 0.6 and 4.9%. With the 
choice N = N2, the corresponding error is located between 12.7 and 42.2%. Results 
for P = 0.4 are presented in table [31 

For p = 0, one obtains b = 1.52, c = 1.26, and d = 1.09. From (fl"32|) with A = 1, 
one obtains b — 1.52, c = 1.23, and d = 1. These values are close to each other as 
expected. 

10. Summary of the results 



The auxiliary field method, which is strongly connected with the envelope theory Pj6], 
is a powerful tool to compute approximate analytical solutions of the Schrodinger 
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Table 3. Eigenvalues e(n, I) of the Hamiltonian (1 1 36 I t with /3 = 0.4, for some 
sets (n,l). First line: value from numerical integration; second line: approximate 
result with N defined by JT37J and JT38J; third line: approximate result with 
N = 2n + l + 3/2. 



1 


e(0,0 


e (i,0 


e(2,J) 


e(3,0 





2.7821 


4.3709 


5.5874 


6.5938 




2.7804 


4.4196 


5.5977 


6.5678 




3.2249 


5.1381 


6.5115 


7.6420 


1 


3.9944 


5.2365 


6.2744 


7.1772 




3.9737 


5.2529 


6.2765 


7.1552 




4.2895 


5.8652 


7.0993 


8.1486 


2 


4.8993 


5.9549 


6.8772 


7.7028 




4.8837 


5.9710 


6.8887 


7.6978 




5.1381 


6.5115 


7.6420 


8.6255 


3 


5.6588 


6.5927 


7.4311 


8.1957 




5.6489 


6.6115 


7.4508 


8.2046 




5.8652 


7.0993 


8.1486 


9.0774 



equation [3J [3]. This method was already used to deal with semirelativistic 
Hamiltonians in previous works [X 91 120[ 126] . In the present paper, we extend this 
technique and apply it to various potentials presenting an interest in atomic and 
hadronic physics: power-law interactions (with special focus on quadratic, linear, and 
Coulomb potentials), square root potential, funnel potential, and Yukawa interaction. 
Both nonrelativistic and ultrarelativistic limits are analyzed. Previous results are 
recovered [12] HH HH] , but new ones are presented. In particular, closed formulas are 
computed for the first time. It is shown that their accuracy can be largely improved by 
slight modifications given by a comparison with exact results coming from a numerical 
analysis. 

It is worth saying that in hadronic and atomic physics, one is sometimes led 
to use running masses, that is basically to deal with a Hamiltonian of the form 
H = (J\/p 2 + m 2 (r) + gU(r), where the mass is position-dependent [39]. Provided 
that an analytical solution, e(m, a), of the Hamiltonian p 2 jm + aU(r) is known, one 
can use the AFM with P(r) = U(r) and V(r) — m 2 (r). The energy formula reads 
vE{v,p) = a 2 + m 2 (I(p)) — pU(I(p)) + ue(u,g + p/v), and analytical solutions can 
then be hoped after elimination of the auxiliary fields for some particular forms of 
m{r). We hope to present explicit applications of this result in future works. 
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Appendix A. Unequal masses 

It is possible that some particular problems require to deal with a system of two 
particles with unequal masses. In this case, a general spinless Salpctcr Hamiltonian 
is given in the rest frame by 




(A.l) 
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A general but useful result can be obtained for Hamiltonians with two different 
masses. Let us consider the following two-body Hamiltonians H = T\ + T% + V, 
Hi = 2Ti + V, and H 2 = 2T 2 + V whose ground state energies are respectively 
E = (<p\H\<f>), E x = {(PxlHilfc), and E 2 = (fa\H 2 \fa). Since H = (H 1 +H 2 )/2, we 
can write 

(<j>\H\<j>) = \((<t>\H 1 \4>) + (ct>\H 2 \4>)). (A.2) 

The Ritz theorem implies that 

E>^(E 1 +E 2 ). (A.3) 

For particular cases, this approximation can be quite good. In [38j . it is shown that 
E ps (Ei + E 2 )/2 for a relativistic Hamiltonian of kind (| A. 1[) with V(r) — ar and 
rrii -C y^a. 

The square roots appearing in the kinetic terms can still be avoided by resorting 
to the AFM as exposed in section [3l But this time, two auxiliary fields, v\ and v 2 , 
have to be introduced. One is led to the Hamiltonian 

H( Vll v 2 ) = + -± + -± + —f r + V(r), A.4 

2 2i>i 2v 2 2M(v\ 1 v 2 ) 

with 

%^)~ (A.5) 
V\ + v 2 

playing the role of a reduced mass. Analytical energy formulas can then be found 
from (|A.4|I , but the minimization of the energy with respect to the two auxiliary fields 
leads to very complicated equations. That is why we have chosen to restrict ourselves 
to the case mi = m 2 = m, implying v\ = v 2 = v = 2M by symmetry. 

A very interesting solvable case corresponds to the situation where the particles 
interact via a Coulomb potential and only one of them is massless. Thus, we are 
concerned with the eigenvalues of the spinless Salpeter equation resulting from the 
Hamiltonian 

a 



H = vT + Vp 2 + 'm 2 • (A.6) 

r 

The introduction of two auxiliary fields fi = v\ and v = v 2 in the nonrelativistic 
Hamiltonian leads to 

2 2v 2M(/j,,v) r 

The original Hamiltonian (|A.6|) is recovered by a proper elimination of /i and v. The 
eigenenergies of (|A.7|) are given by 

v) = -^p M ^> + ^ + ^ ( A -8) 

with N = iV_ i = n + I + 1 a priori. 

It is easier to begin the minimization by calculating the derivative with respect to 
fj,, dE/d^i — 0, using the property OM/dfJ, = M 2 /fi 2 . One obtains the minimal value 

Mo(f) = (jf - l) v. (A.9) 

Reporting in E(jjl (v),v) = v[a/N - a 2 j(2N 2 )\ + m 2 j(2v), one is led to the 
minimization of this function of v only. It is easy to find the value vq which minimizes 
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E, and the value of the energy E — _E(/io(fo), vq) which corresponds to this minimum. 
Explicitly, one obtains 



E = 2mJ—(l-—). (A.10) 
V 2N V 2NJ V ' 

In this formula, one has in principle N = n + I + 1, but, as we pointed out several 

times, it is justified to take a more sophisticated expression in order to get a better 

accuracy. It can be seen from (IA.10[) that rim m _,o E = 0, as expected. 

Appendix B. Some polynomial equations 

Finding analytical energy formulas for the potentials that we study in this work 
requires an analytical knowledge of the roots of particular cubic and quartic equations. 
We sum up these equations in this appendix and put their roots in a form that is as 
convenient as possible to deal with. Notice that the needed polynomial equations are 
the same as those appearing in our previous work [3]. 
We begin by the cubic equation (Y > 0) 

x 3 ± 3a; - 2Y = 0, (B.l) 

for which there exists only one positive root given by 

1/3 / \-i/3 



F±(Y) = [Y + y/Y 2 ±l) T [Y+y/Y 2 ±l) . (B.2) 



Written in the above form, it seems that F-(Y) is not properly defined for Y < 1. 
But, for this range of Y values, one can show that 

F_ (Y) = 2 cos Q arccos Y^j . (B.3) 

So F—(Y) is well defined for all positive values of its argument. It can be checked that 
the following approximate forms hold 

2Y 

F + (Y)^— if y«l, (B.4) 

F_{Y)^V3 + \ if (B.5) 

F±(Y) « (2Y) 1/3 if F>1. (B.6) 
The quartic equation which gives the most pleasant form for the roots is (Y > 0) 
4a; 4 ± 8a; - 3F = 0. (B.7) 
For each sign, there exists only one positive root given by 

G±(Y) = T ^^v(Y) + y^V(Y))-^ V(Y), (B.8) 

with 

1/3 / \-i/3 



V(Y) = (2 + V4 + F3) -Y (2 + V4 + F3) . (B.9) 



The following approximate expressions can also be used for simplicity 

3F 

G+(r)«— if (b.io) 

G_(Y)«2 1/3 + — if F<1, (B.ll) 
8 

/oy\ 1/4 

G±(Y) « i- if Y»l. (B.12) 
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